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ABSTRACT 

Aims. In 2010, detailed observations have been published that seem to indicate another super-Earth planet in the 
system of Gliese 581 located in the midst of the stellar climatological habitable zone. The mass of the planet, known 
as Gl 581g, has been estimated to be between 3.1 and 4.3 Af®. In this study, we investigate the habitability of Gl 581g 
based on a previously used concept that explores its long-term possibility of photosynthetic biomass production, which 
has already been used to gauge the principal possibility of life regarding the super-Earths Gl 581c and Gl 581d. 
Methods. A thermal evolution model for super-Earths is used to calculate the sources and sinks of atmospheric carbon 
dioxide. The habitable zone is determined by the limits of photosynthetic biological productivity on the planetary 
surface. Models with different ratios of land / ocean coverage are pursued. 

Results. The maximum time span for habitable conditions is attained for water worlds at a position of about 0.14 ± 
0.015 AU, which deviates by just a few percent (depending on the adopted stellar luminosity) from the actual position 
of Gl 581g, an estimate that does however not reflect systematic uncertainties inherent in our model. Therefore, in 
the framework of our model an almost perfect Goldilock position is realized. The existence of habitability is found to 
critically depend on the relative planetary continental area, lending a considerable advantage to the possibility of life 
if Gl 581g's ocean coverage is relatively high. 

Conclusions. Our results are a further step toward identifying the possibility of life beyond the Solar System, especially 
concerning super-Earth planets, which appear to be more abundant than previously surmised. 
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lar climatological habitable zone. Its mass is estimated to 

The nearby M-type dwarf Gliese 581 is an important tar- be between 3>1 and 4.3 Mffi) noting that the p uta tive upper 

get of previous and ongoing planetary search efforts, and mass limit is partially constrained by the assumption that 

so far four, but possibly up to six planets around this star none of the p i ane tary orbital eccentricities exceeds 0.2. The 

have been discovered. It is noteworthy that until recently distance of Gl 581g is estimated to be 0. 14601±0. 00014 AU, 

the possible existence of Earth-type and super-Earth plan- and its eccentrici ty is almost ide n tical t o zero, 

ets in extrasplar systems has beenjiigjriy speculative (e.g., According to [Valencia "et al.1 (l20f)fih . super-Earths are 

i Butler et al. || 2006| ).| Bonfils et al. | (|2005|) reported the detec- rQcky planets from Qne tQ about ten Earth masseg with a 

■ tion of a Neptune-size planet around Gl 581 (M3 V) . Gl 581 chemical and mineral composition akin to that of Earth. In 

r ^ , is at a distance of 6.26 pc; it has a mass of about 0.31 M Q the followingj we adopt the hypothesis that this is indeed 

_Cd_, and a lumin osity of 0.01 3 ± 0.002 L e (see below). the cas6; and considcr for G1 581g our mode l previously 

Later on,|Udry et aTJ (|2007|) announced the detection of developed for the Earth, using scaling laws if appropriate, 

two so-called "super-Earth" planets in this system, Gl 581c T hi s ap proach has already b een applied to Gl 581c and 

and Gl 581d, with minimum masses of 5.06 and 8.3 M e and Ql 581d (|von Bloh et al.ll2007t) . as well as to simulations of 

with semi-major axes of 0.073 and 0.25 AU, respectively. In hypothetical super-Earth planets in various other systems 

2009, Gl 581e was discovered, another super-Earth planet i.e., 47 UMa and 55 Cnc (jCuntz et aT1l2003l: iFranck et all 

([Mayor et alj|2009|). Meanwhile the position of Gl 581d has [200I Ivon Bloh et"aT] [2003|) as well as systems u ndergoing 

been r e-identified as 0.22 AU. Finally, in 2010, |Vogt et al.| re d giant branch evolution (jvon Bloh et al.|[2009l ). 

42Q10D reported a still controversial discovery. The authors xhe main question ig whetner G1 581gj if existing, of- 

argue in favor of the existence of two additional super-Earth fcrs thc p rinc i pa i possibility of life, i.e., whether it lies 

planets, namely Gl 581f and Gl 581g. Note that Gl 581g is wit hin the stellar habitable zone (HZ), as already argued by 

of particular interest to thc study of planetary science and |y ogt et aL | (goH|. Typically, stellar HZs are defined as re- 

gions around the central star where the physical conditions 

Send offprint requests to: W. von Bloh are favourable for liquid water to exist at the planet's sur- 
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face fo r a period of time suff icient for biological evolution to 
occur. iKasting et al.l (|l993j) calculated the HZ boundaries 
for the luminosity and effective temperature of the present 
Sun as R iu = 0.82 AU and R out = 1.62 AU. They defined 
the HZ of an Earth-like planet as the region where liquid 
water is present at the surface. 

According to this definition, the inner boundary of the 
HZ is determined by the loss of water via photolysis and hy- 
drogen escape. The outer boundary of the HZ is determined 
by the condensation of CO2 crystals out of the atmosphere 
that attenuate the incident sunlight by Mie scattering. The 
critical CO2 partial pressure for the onset of this effect is 
about 5 to 6 bar. Howeve r, the cooling effect of CO 2 cloud s 
has been challenged by iForget fe Pierrehumbertl (|1997|) . 
CO2 clouds have the additional effect of reflecting the out- 
going thermal radiation back to the surface. The precise in- 
ner and outer limits of the climatic habitable zone are still 
unknown due to the limitat i ons of the existing climate mod- 
els. Recently, iHeller et al.l (|2011h extended the concept of 
the HZ by including constraints arising from tidal processes 
due to the planet's spin orientation and rate. Effects caused 
by tilt erosion, tidal heating, and tidal equilibrium rotation 
seem to be especially important for potentially habitable 
planets (as Gl 581d and Gl 581g) around low-mass stars. 



For limit ations of the planeta r y habitability of M-type stars 
see, e.g..lTarter et al.l (|2007t ). iGuinan fc Engld (|2009D . and 



ISegura et all (|201lh . 

The luminosity and age of the central star play im- 
portant roles in the manifestation of habitability. The lu- 
minosity of Gl 581 can be obtaine d by (1) photometry 
(jBonfils et al.ll2005l: lUdrv et al.ll2007[ ) and (2) th e applica- 
tion of the mass-radius relationship (|Ribasll20"06t ) together 
with the spectroscopicall y determined stell ar effective tem- 
perature of T e = 3480 K dBean et al.ll20q6[). Bot h methods 
yield L = 0.013 ± 0.002 L ? . ISelsis et al.l (|2007l) estimated 
the stellar age as at least 7 Gyr based on the non-detection 
of Gl 58 1's X-ray flux considering the sensitiv ity limit of 
ROSAT ilSclm.it t et al.lll99fil IVoges et al.ll2000ri . 

In the follow ing, we adopt a defin ition of the HZ pre- 
viously used by iFranck et al.l (|2000al fbh . Here habitability 
at all times does not just depend on the parameters of the 
central star, but also on the properties of the planet. In 
particular, habitability is linked to the photosynthetic ac- 
tivity of the planet, which in turn depends on the planetary 
atmospheric CO2 concentration together with the presence 
of liquid water, and is thus strongly influenced by the plan- 
etary dynamics. We call this definition the photosynthesis- 
sustaining habitable zone, pHZ. In principle, this leads to 
additional spatial and temporal limitations of habitability 
because the pHZ (defined for a specific type of planet) be- 
comes narrower with time owing to the persistent decrease 
of the planetary atmospheric CO2 concentration. 



2. Habitability of super-Earth planets 

2.1. Photosynthesis-sustaining habitable zone (pHZ) 

The climatic habitable zone at a given time for a star with 
luminosity L and effective temper ature T e different from the 
Sun can be calculated foll owing Underwood et all (|2003l ) 
based on previous work by IKasting et al.l (|1993l ) as 
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Fig. 1. Earth system box model. The arrows indicate the 
different forcing and feedback mechanisms. The bold ar- 
rows indicate negative feedback operating towards climate 
stabilization. 



with Si n (T e ) and S ou t(T e ) described as second order poly- 
nomials. 

To assess the habitability of a terrestrial planet, an 
Earth-system model is applied to calculate the evolution of 
the temperature and atmospheric CO2 concentration. On 
Earth, the carbonate-silicate cycle is the crucial element for 
a long-term homeostasis under increasing solar luminosity. 
On geological time-scales, the deeper parts of the Earth are 
considerable sinks and sources of carbon. 

Our numerical model previ ously applied to Gl 581c and 
Gl 58 Id (jvon Bloh et al.|[2~007|) couples the stellar luminos- 
ity L, the silicate-rock weathering rate F wr and the global 
energy balance to obtain estimates of the partial pressure of 
atmospheric carbon dioxide Pco 2 1 the mean global surface 
temperature T sur f, and the biological productivity n as a 
function of time t (Fig.QJ. The main point is the persistent 
balance between the CO2 sink in the atmosphere-ocean 
system and the metamorphic (plate tectonic) sources. This 
is expressed through the dimensionless quantities 



U(t) ■ f A (t) = /„(*), 



(2) 



where / wr (t) = F m (t)/F wr fi is the weathering rate, /a(^) = 
A c (t)/A C fl is the continental area, and f sr (t) = S(t)/So 
is the areal spreading rate, which are all normalized by 
their present values of Earth. Eq. ^ can be re arranged b 
introducing the geophysical forcing ratio GFR (|Volklll98 
as 



/wr(T E 



surf, PCQ 2 ) 



fsr 

Ta 



GFR(t). 



(3) 



Here we assume that the weathering rate only depends on 
the global surface temperature and the atmospheric CO2 
concentration. For the investigation of a super-Earth under 
external forcing, we adopt a model planet with a prescribed 
continental area. The fraction of continental area with re- 
spect to the total planetary surface /a is varied between 
0.1 and 0.9. 
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Table 1. Size and gravity parameters of the Gl 581g models 
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Parameter 




Value 




Unit 


Description 






3.1 M e 


4.3M® 






a 


9.81 


16.5 


19.2 


m s - ^ 


gravitational acceleration 


Rp 


6378 


8657 


9456 


m 


planetary radius 


R c 


3471 


4711 


5146 


in 


inner radius of the mantle 


Rm 


6271 


8511 


9298 


m 


outer radius of the mantle 



The connection between the stellar parameters and the 
planetary clim ate can be obtai ned by using a radiation bal- 
ance equation (|Williamslll998h 



L 



4ttR 2 



[1 - a(T surf ,Pco 2 )] = 4/ fl (T sur f,P C o 2 ), 



(4) 



where a denotes the planetary albedo, Ir the outgoing in- 
frared flux, and R the distance from the central star. The 
climate model does not include clouds, which are particu- 
larly hnp_OTtairt for_determining the inner boundary of the 
HZ (jSelsis et al J 120071 ). The Eqs. © and (0} constitute a 
set of two coupled equations with two unknowns, T sur f and 
Pco 2 i if t ne parameterization of the weathering rate, the 
luminosity, the distance to the central star and the geo- 
physical forcing ratio are specified. Therefore, a numerical 
solution can be attained in a straightforward manner. 

The photosynthesis-sustaining HZ around Gl 581 is de- 
fined as the spatial domain of all distances R from the cen- 
tral star where the biological productivity is greater than 
zero, i.e., 



pHZ := {R\IL(Pco 2 (R,t),T snrl (R,t))>0}. 



(5) 



In our model, biological productivity is considered to be 
solely a function of the surface temperature and the CO2 
partial pressure in the atmosphere. Our parameterization 
yields zero productivi ty for T surf £ 0°C o r T surf > 100°C 
or P C o 2 < 1CT 5 bar (|Franck et al.ll2000al) . The inner and 
outer boundaries of the pHZ do not depend on the detailed 
parameterization of the biological productivity within the 
temperature and pressure tolerance window. 

2.2. Comments on the thermal evolution model 

Parameterized convection models are the simplest models 
for investigating the thermal evolution of terrestrial planets 
and satellites. They have been successfully applied to the 
evolution of Mer cury, Venus, Earth, Mars, and the Moon 
([Stevenson et all 119831: ISleed |2000D . iFranck fc Bounamal 
(|l995l) studied the thermal and volatile history of Earth and 
Venus in the framework of comparative planetology. The in- 
ternal structure of massive terrestrial plan ets with one to 
ten E arth masses has been investigated by IValencia et al.l 
(|2006l) to obtain scaling laws for total radius, mantle thick- 
ness, core size, and average density as a function of mass. 
Similar scaling laws were found for different compositions. 
We will use these scaling laws for the mass-dependent prop- 
erties of super-Earths an d also the mass-independen t ma- 
terial properties given by IFranck fc Bounamal (jl995l ) . 

The thermal history and future of a super-Earth has 
to be determined to calculate the spreading rate for solving 
key Eq. @. A parameterized model of whole mantle convec- 
tion including the volatile exchange between the mantle and 
surface reservoirs (jFranck fc Bounamal [l995t lFranckj|1998[ ) 
is applied. The key equations used in our present study are 



in accord wi th our previous work focused on Gl 581c and 
Gl 581d; see Ivon Bloh et al.l (|2007l ). for details. A key ele- 
ment is the computation of the areal spreading rate S; note 
that S is a function of the average mantle temperature T m , 
the surface temperature T sur f , the hea t flow from the mantle 
(7 m , a nd the area of ocean basins Aq (|Turcotte fc Schubert! 
120021) . It is given as 



.5' 



q^wKAo 



4 k ^ ( T m ur f ) 2 



(0) 



where k is the thermal diffusivity and k the thermal con- 
ductivity. To calculate the spreading rate, the thermal evo- 
lution of the mantle has be to computed: 



-7rpc(R 3 m -R 3 c ) 



dT„ 



-4nRl l q ri 



+ -7rE(t)(R 3 m -R 3 c ),(7) 



3 r-,- m ~c, dt -—to,.,. ■ 3 

where p is the density, c is the specific heat at constant 
pressure, E is the energy production rate by decay of ra- 
diogenic heat sources in the mantle per unit volume, and 
R m and R c are the outer and inner radii of the mantle, re- 
spectively. To calculate the thermal evolution for a planet 
with several Earth masses, i.e., 3.1 and 4.3 M$ as pursued 
for Gl 581g in our present study, the planetary parameters 
have to be adjusted. Thus we assume 

0.27 

(8) 



R 



M 

Mm 



where R p is the planetary radius, see I Valencia et afl (|2006l ). 
The total radius, mantle thickness, core size, and average 
density are all functions of mass, with the subscript © de- 
noting Earth values. The exponent of 0.27 has been ob- 
tained for super-Earths. The values of R p , R m , R c , as well 
as the other planetary properties are scaled accordingly. It 
means that R p , R m and R c increase by a factor of 1.36 for 
M = 3.1M e and 1.48 for M = 4.3M®. 

Table [T] gives a summary of the size parameters for the 
models of the p l anets with 3.1 and 4.3 M e ; see the study by 
Ivon Bloh et al.l (|2007|) for additional information. The val- 
ues for an Earth-size planet are included for comparison. 
The onset of plate tectonics on massive terrestrial planets 
is a topic of controversy. While [O'Neill fc Lenardid (|2007l) 
stated that there might be in an ep isodic or stagnant lid 
regime, IValencia fc O'Connelll (|2009l) proposed that a more 
massive planet is likely to convect in a plate tectonic regime 
similar to Earth. Thus, the more massive the planet is, 
the higher the Rayleigh number that controls convection, 
the thinner the top boundary layer (lithosphcre), and the 
higher the convective velocities. In the framework of our 
model, a platc-tcctonic-driven carbon cycle is considered 
necessary for carb on-based life. This app roach follows the 
previous work by Ivon Bloh et al.l (|2007t ). who gave a de- 
tailed discussion of the equations and parameters used in 
their study for the super-Earths with 5 and 8 M®, identified 
as Gl 581c and Gl 581d, respectively. 
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Fig. 2. The pHZ of Gl 581 for a super-Earth of (a) M = 
3.1 A/© and (b) M = 4.3M© as a function of planetary 
age. The relative continental area is varied from 0.1 to 0.9. 
The stellar l umino sity is assumed to be O.O13L0. Following 
ISelsis et "ail (|2007l ) and references therein, the stellar age is 
at least 7 Gyr, indicated by a dashed line. For comparison, 
we show the positions of Venus, Earth, and Mars scaled to 
the luminosity of Gl 581. Note that the position of Gl 581g 
almost exactly coincides with the luminosity-scaled position 
of Earth. 



3. Results and discussion 

The photosynthesis-sustaining habitable zone is calculated 
for super-Earth planets with 3.1 and 4.3 M© and the re- 
sults are depicted in Fig. [2] The maximum life span of a 
biosphere is given at the point in time when the pHZ van- 
ishes. This maximum life span strongly depends on the rel- 
ative continental area r and increases with decreasing r. 
Therefore, "water worlds" are favored in the facilitation of 
habitability as previously obtained in model s of fictitious 
super-Earth planets for 47 UMa and 55 Cnc (jFranck et al.l 
l2003t Ivon Bloh et al.ll2003l ). In this context, water worlds 
are planets of non-vanishing continental area mostly cov- 
ered by oceans. The climate of a planet fully covered by 
oceans is not stabilized by the carbonate-silicate cycle. The 
maximum life span for a water world with r = 0.1 and a 
planetary mass of M = 3.1M® is 15.4 Gyr, whereas it is 
16.3 Gyr for M = 4.3M©. Both times are considerably 
longer than the estimat ed age limit of Gl 581, which is 
7 Gyr (ISelsis et alJl2007h . 
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Fig. 3. Life span for Gl 581g with M = 3.1M® (+) and 
M = 4.3M e ( function of the relative continental 

area. For comparison the dashed curve denotes the life span 
of an Earth- mass planet. The horizontal line indicates the 
lower limit for the age estimate of 7 Gyr. 

This maximum life span can be realized for a super- 
Earth at a distance from the central star of about 0.14 ± 
0.015 AU. This uncertainty is reflective of the uncertainty 
in the luminosity estimate of Gl 581, noting however that 
our models are also subject to various systematic uncertain- 
ties (see Sect. 2). A full assessment of those uncertainties, 
which are expected to be larger than the impact of the un- 
certainty in the stellar luminosity, can best be obtained by a 
comparison with future alternative models. It is noteworthy 
that the orbital position of Gl 581g at 0.146 AU, if existing, 
is largely consistent with this ideal position. Thus, we can 
conclude that Gl 58 lg has an almost perfect Goldilock po- 
sition in the star-planet system of Gl 581, and it is the best 
candidate for extra-terrestrial habitability so far. Because 
of the Goldilock position, located well within the HZ, clouds 
play no major role for our results. Any shift (within limits) 
of the inner or outer boundary of the HZ will not affect the 
maximum life span obtained by our model. The optimum 
position is not only the geometric mean of the pHZ, but 
it is also the position where the maximum life span of the 
biosphere is realized. 

But the pHZ also strongly depends on the planetary age. 
Figure [3] depicts the life span for a super-Earth with 3.1 and 
4.3 M© at 0.146 AU as a function of the relative continental 
area r. For a relative continental area larger than 0.6, the 
realized life span is shorter than 7 Gyr, which is the esti- 
mated lower limit for the stellar age of Gl 581 star-planet 
system. In this case, no habitable conditions on Gl 58 lg 
would exist in the framework of the adopted geodynamic 
model. If we place an Earth twin (r ~ 0.29) at the orbital 
position of Gl 581g, habitable conditions would cease just 
at an age of 7 Gyr. In conclusion, we have to await future 
missions to identify the pertinent geodynamical features of 
Gl 581g (in case it does exist) and to search for biosigna- 
tures in its atmosphere to gain insight into whether or not 
Gl 581g harbors life. 
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